Photoinduced coherent oscillations in the one-dimensional two-orbital Hubbard model 
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We study photoinduced ultrafast coherent oscillations originating from orbital degrees of freedom 
in the one-dimensional two-orbital Hubbard model. By solving the time-dependent Schrodinger 
equation for the numerically exact many-electron wave function, we obtain time-dependent optical 
response functions. The calculated spectra show characteristic coherent oscillations that vary with 
the frequency of probe light. A simple analysis for the dominant oscillating components clarifies 
that these photoinduced oscillations are caused by the quantum interference between photogenerated 
states. The oscillation attributed to the Raman-active orbital excitations (orbitons) clearly appears 
around the charge-transfer peak. 

PACS numbers: 78.20.Bh, 71.10.Fd, 75.25. Dk, 78.47.J- 



Photoinduccd phenomena of strongly correlated elec- 
tron systems have attracted much attention recently^— 
For example, there have been many studies on photoin- 
duced macroscopic changes in electronic states, often 
called "photoinduced phase transitions" (PIPTs)4^— 

These photoinduced phenomena often accompany sub- 
sequent non-equilibrium dynamics. One typical example 
is coherent oscillations observed after the rapid photoin- 
duced changesJ^— These oscillations involve much in- 
formation of characteristic collective modes of the sys- 
tems, phonon, orbiton, and so on. Hence investigating 
the coherent oscillations provides us insight into roles of 
these modes in the photoinduced phenomena. 

Until a few years ago, experimental studies have used 
relatively long pulses (~ lOOfs), which allow us to detect 
only slow lattice dynamics However, recent develop- 
ment of experimental technique that provides sub-10-fs 
pulses enables us to observe much faster dynamics. In 
particular, considerable experimental effort has been de- 
voted to the study of the ultrafast oscillations in tran- 
sition metal oxides, which have fast vibrational phonon 
modeSfi^— or orbital excitations^ 

In contrast to these experimental achievements, theo- 
retical studies on the photoinduced ultrafast oscillations 
have not been carried out so intensively^^— Although 
some of the authors and coworkers have provided a theo- 
retical description on dynamics of an organic compound 
(EDO-TTF) 2 PF 6 , the treatment for the lattice degrees of 
freedom is limited to a classical onei 18 ' 19 A quantum the- 
ory for the same material^ where quantized phonons are 
dealt with, focuses on the slow lattice dynamics. Thus al- 
ternative quantum-mechanical treatment is needed to de- 
scribe the ultrafast oscillations of excitations with much 
higher frequencies. 

In this Rapid Communication, we present a theoretical 
study of the photoinduced ultrafast coherent oscillations 
of the one-dimensional (ID) two-orbital Hubbard model 
coupled with static lattice distortion, which is a ID ana- 



log of transition- metal oxides with orbital degrees of free- 
dom. Numerically calculated time-dependent optical re- 
sponse functions show clear ultrafast coherent oscillations 
that vary with the frequency of probe light. An analysis 
of optical excitation processes contributing the dominant 
oscillating components clarifies that (i) the photoinduced 
coherent oscillations are caused by the quantum interfer- 
ence between eigenstates included in the photoexcited 
state and (ii) the oscillation around the charge-transfer 
(CT) peak results from the Raman-active two-orbiton 
state while the oscillation in the low-energy region is 
caused by a onc-holon-doublon (hd)-pair excitation. The 
excitation process for the former case is illustrated in 

Fig.ru 




FIG. 1: The excitation process contributing to the coherent 
oscillation in the CT region. The bold solid arrows show 
the optical transitions connecting the initial antiferro-orbital 
state (1) and the final two-orbiton state (3) via the interme- 
diate optically excited state (2). The quantum interference 
between (1) and (3) generates the coherent oscillation with 
the frequency A20, the energy gap between the two states. A 
clear oscillation appears in the case where the probe frequency 
cjprb is almost equal to the optical gap A opt . 



In this work, we use the ID two-orbital Hubbard model 
coupled with static lattice distortion. The Hamiltonian 
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is given by 
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the transient optical response functio: 
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t, |) at orbital r(= 1,2) at site 
niri, and Qz is the Jahn- 



wherc cj Ta (ci Ta ) is the creation (annihilation) operator of 
an electron with spin a 
I, n iTa = c\ Ta c lTa , m T = ni T ^ 
Teller- type lattice distortion. U, U', J, and J' denote 
intraorbital Coulomb, interorbital Coulomb, intcrorbital 
spin exchange, and intcrorbital pair hopping interactions, 
respectively. We also note that the following relations 
U = U' + 2 J and J' = J hold M- The electron-lattice 
coupling and the elastic constant are given by g and K, 
respectively We treat the quarter-filled JV-site chain with 
N = 4 and impose the periodic boundary condition. 

The timc(r) -dependent transfer integral t(r), which is 
finite only between the same orbitals of neighboring sites, 
is introduced as i(r) = toe t ^ ae ' ,ic ' A ^ T > , where to is the bare 
transfer integral, e is the absolute value of the electronic 
charge, a is the lattice spacing, and c is the velocity of 
light. In the following, we use the unit to = e = a = c = 
h = 1. The pump laser pulse is represented by the vector 
potential A(t) given by 
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where F is the amplitude of the electric field, r c is the 
central time of the pump field, and T pmp defines the width 
of the Gaussian function. We set t c = 10 and T pmp = 1. 
When we set the bare transfer integral to = 0.1 eV, the 
pulse width 2T pmp corresponds to about 13 fs, which is 
same order of the pulse width of recent experiments 
The frequency w pmp is set to the optical gap A opt . 

The procedure of calculation is as follows. First of all, 
we obtain the ground state |<^o) and the stable lattice 
distortion Qi with no pump field by iterative applica- 
tion of the Lanczos diagonalization and the Hellmann- 
Feynmann theorem^ d{(f>o\H\(f>o) /dQi = for all I. The 
obtained stable configuration is found to be the staggered 
distortion Qi = {-l) l Q st for I = [0, • • ■ ,N-1]. Then, we 
calculate the state \4>(r)) by solving the time-dependent 
Schrodinger equation i4-\^) = V) • The ground state 
\4>o) is used as the initial state and the lattice configu- 
ration is fixed. The Schrodinger equation is numerically 
solved by expanding the exponential time evolution op- 
erator with time slice dr=0.02J^ 

To observe the time-dependent dynamics, we calculate 
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given by 



-Mr)), 



iS + E-Ho 

(3) 

where j = it J2i T a(4Ta c Wra-cj +lTa ci Ta ) is the current 
operator^ S is a broadening parameter set at 2.0, Ho is 
Hamiltonian (JTJ with t(r) = t , and E = (f/j(r)|'Ho|V'W}- 
We also note that A opt is obtained from the lowest peak 
of J(o; P rb, t) at t = 0. 

We set other parameters as U' = 20, J = 5, g = 0.4, 
and K = 1. For these parameters the ground state is 
in the ferromagnetic phase with finite lattice distortion 
Q s t ~ 0.36. Because we focus on the photoinduced coher- 
ent oscillation phenomena in this work, the pump field is 
set weak, F=2.0, which does not cause a PIPT. 
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FIG. 2: (Color online). The optical response function 
iCwprbjT) for (a) Ljp r b = 10 and 20, and (b) for oj plh = 1. 
The inset of (b) shows the optical response I(w pI b, r = 0). 

Figure^ (a) shows the optical response I(w pr b,r) for 
Wp r b = 10 and 20, which are the lower and higher sides of 
the CT peak [see the inset of Fig. [2] (b)]. After the drastic 
change induced by the pump field, I(a> pr b = 10,20, r) 
shows clear coherent oscillations. Their periods are about 
8 for both w pr b=10 and 20, which suggests that their 
origins are the same low-lying excitations. By contrast, 
the oscillation for I(w pr b = 1,t) with period ~ 2 [see 
Fig. [2] (b)] is evidently caused by other excitations. 

To clarify the origins of the coherent oscillations 
around the CT peak, we calculate the Fourier transform 
of 7(a; pr b,T) and two spectral functions defined below. 
The Fourier transform in the time domain re [t,, r e ] is 
given by 
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One of the spectral functions is that detects Raman- 
active excitations, defined by 
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where \jj) = jj\<f> ) - |0 o )(0o|Bl0o)- The other one is 
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the orbital dynamical structure factor T z (q,ui) given by 



TV H u> + 10 + e - no 

where T z = £ Tfe~ iql and Tf = (n n - n /2 )/2. 
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FIG. 3: (Color online). The Fourier transforms 7(w P rb,w) for 
^ P rb = 10 and 20 in the time domain [Ti,T e ] = [20, 100], and 
the spectral functions x{uj) and T z (q = n, u) . The broadening 
parameter 8' is set at 0.05. 

Figure El shows I(u> pT b, uj) for the time domain [n, r e ] = 
[20, 100], and T z (q = tt,u>). All the functions have 

two distinct peaks: the dominant one with frequency ui = 
0.8, and the sub-dominant one with ui = 1.4. Hence we 
conclude that the two peaks of I(u; pro bc, w) correspond 
to the Raman-active orbital excitations. These orbital 
excitations are described by the effective model for the 
orbital degrees of freedom 2 ^ given by 



U eS = JJ^fr f l+1 - H cff ^(-l)'Tf 
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where J is the antiferro-orbital superexchange constant, 
J = 4fo/(£/' ^ J), and H c g is the effective staggered 
field, -ffeff = 2(7<3st- The spin part is omitted since the 
system is ferromagnetic. For the two-orbital system, the 
Raman-active orbital excitations have been discussed and 
shown that they are two-orbiton processes with a finite 
excitation gap due to the staggered fields 
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than the frequencies of the orbital excitations. The ori- 
gin is elucidated by calculating an analog of x(w) for the 
lowest optical excitation |c/> opt ) defined by 
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where \jj') = jj\4> op t) ~ |0o P t)(0o P tbj|0o P t) and e opt is 
the energy of 1 O pt) ■ The state |0 O pt) i s the lowest one- 
holon-doublon(hd)-pair excitation ) 27 ' 28 which gives the 
main contribution of the low-w pr b component after pho- 
toexcitation. x( w )o P t shown in Fig. |4] has a clear peak 
at ui = 3.1, which is at the same location of that of 
-f(w pr b = The state corresponding to this peak 

is another one-hd-pair excitation | lhd') in Fig. [5] (for de- 
tails, see below), taking account of its eigenenergy; it is 
higher than that of |<^ pt) by the order of to and much 
lower than the two-hd-pair excitations. 

Now, let us discuss the reason why the different oscil- 
lations are observed by changing w pr b- To this end, we 
expand the quantum state \iP(t)) as follows; 
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where \a) is an eigenstate of Ho and has the eigenenergy 
e a . We now assume that the total energy E is almost 
equal to the ground-state energy eo since the pump field is 
weak in this work, and we thereby obtain the expression 
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x (7b1/3)(/S|j|«)/ L Krb + eo-^), (10) 

where fh(x) is the Lorcntzian function = x 2^_g2 ■ 

Equation fp~0|) tells us the following points: (i) there are 
three important states, the initial state the final state 
I7), and the virtually excited state |/3), which are con- 
nected by the matrix element of j. (ii) The coherent 
oscillation occurs as a quantum interference between \a) 
and 1 7) and its frequency is equal to the energy difference 
e 7 — e a . (iii) The oscillation appears for w pr b ~ ep — eo- 
By using these points, we discuss the coherent oscilla- 
tion around the CT gap, i.e., w pr b ~ A opt = e opt — eo- In 
this case, the relevant virtual state \/3) is \4> pt), and the 
important initial state is the ground state \4>o)- Then the 
final state is expected to be the two-orbiton state |20), 
which is detected as the main peak of T z (q = n,u>). The 
schematic picture of this transition process is shown in 
Fig. [5] (a), and the dominant component is given by 



FIG. 4: (Color online). The Fourier transform I(u} pl h,u}) for 
i^prb = 1 in the time domain [Ti,T e ] = [20, 100] and x( w )o P t 
with 8' = 0.05. 

Now we turn our attention to the low-w pr b region. The 
Fourier transform for uj pr h = 1, depicted in Fig. 21 dis- 
plays the dominant peak at ui = 3.1, which is much higher 



/(wprb.r) ~ ^C; o C e lA ^ T (2O\M opt ) 

X (0o P t|j|0O>/L(w prb - A opt ) + c.c. ,(11) 

where A 2 o = £20 ~ e o- This expression clearly shows that 
the two-orbiton excitation |20) is observed as the main 
oscillating component and that |20) is Raman-active. 
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Now, let us discuss the low-Wprb region. In this region 
the relevant initial state is the lowest optical excitation 
|0opt)- Then the main virtual state is the two-orbiton 
state 1 20), and the final state should be another one-hd- 
pair state | lhd') [sec Fig [5] (b)]. We thereby obtain the 
main contribution; 

x (2O|j|0 op t}/ L (wprb - A 20 ) + c.c. . (12) 



From Eq. (|12j) . we can see that the frequency shown in 
Fig. @] is equal to the gap between the lowest one-hd- 
pair state |</> pt) and another one-hd-pair state |lhd). 
In general, one-hd-pair states form a continuum in the 
thermodynamic limit ^ implying that an infinite number 
of coherent oscillations can contribute to /(oj pr b,T). As 
a result, the coherent oscillation in the low-Wprb region 
may disappear because of the superposition of the infinite 
oscillating components. 
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FIG. 5: The excitation processes contributing to the coherent 
oscillations (a) in the CT region and (b) in the low-w pr b re- 
gion. The solid arrows show the optical transition processes 
from the initial state |a) to the virtual state |/3) and the 
dashed arrows show those from \/3) to the final state I7). 

Here we discuss the experimental realization of the 
photoinduced oscillations caused by the orbital excita- 
tions. The appearance of the photoinduced oscillation 
caused by the orbital excitations is suggested on the basis 
of the experimental results for Pro.7Cao.3Mn03p^ a Mn 



perovskite with three-dimensional structure. However, 
there is a puzzling fact that the oscillation is observed 
only above the orbital melting temperature. Here, we 
note that our theory does not prohibit such an oscillation 
if Raman-active orbital excitations exist in the disordered 
phase. In addition, if other excitations such as phonons 
have dominant Raman intensity, it might be difficult to 
distinguish the orbital excitations by using the Fourier 
transformation even in the ordered phase. 

A more appropriate candidate of the quasi- ID system 
is LaVO3 ) 30 i 31 where the Raman- active two-orbiton ex- 
citations exist^ As for the photoinduced properties, no 
coherent oscillation has been detected while a photoin- 
duced Drude-type spectral weight has been observed??-^ 
Nonetheless, the femtosecond time-resolved reflection 
spectroscopy would clarify the oscillations because there 
found the Raman-active orbital excitations with the fre- 
quencies of 43 and 62 meV^ which correspond to os- 
cillations with the time-periods of 96 and 67 fs. Other 
quasi-lD materials, including KCuFa^ 3 - would be alter- 
native candidates with orbital degrees of freedom. 

In summary, our quantum-mechanical treatment pro- 
vides a simple picture for the photoinduced ultrafast co- 
herent oscillations; the oscillations observed in the optical 
response are caused by the quantum interference between 
the eigenstates included in the photoexcited state. The 
difference of the virtual optical excitation process results 
in the o> pr b dependence of the oscillations. 
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